Solution of a statistical mechanics model for pulse formation in lasers 



o 
o 

(N 

> 
O 



Omri Gat0 Ariel GordonQ and Baruch Fischcjj] 
Department of Electrical Engineering, Technion, Haifa 32000, Israel 
(Dated: February 6, 2008) 

We present a rigorous statistical-mechanics theory of nonlinear many mode laser systems. An im- 
portant example is the passively mode-locked laser that promotes pulse operation when a saturable 
absorber is placed in the cavity. It was shown by Gordon and Fischer [l| that pulse formation is 
a first-order phase transition of spontaneous ordering of modes in an effective "thermodynamic" 
system, in which intracavity noise level is the effective temperature. In this paper we present a 
rigorous solution of a model of passive mode locking. We show that the thermodynamics depends 
on a single parameter, and calculate exactly the mode-locking point. We find the phase diagram 
and calculate statistical quantities, including the dependence of the intracavity power on the gain 
saturation function, and finite size corrections near the transition point. We show that the thermo- 
dynamics is independent of the gain saturation mechanism and that it is correctly reproduced by a 
mean field calculation. The outcome is a new solvable statistical mechanics system with an unstable 
self-interaction accompanied by a natural global power constraint, and an exact description of an 
important many mode laser system. 

PACS numbers: 42.55.Ah, 42.65.-k, 05.70.Fh 
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I. INTRODUCTION 

Lasers can produce light in continuous wave (cw) or 
pulsed manners. A special pulsed operation is mode- 
locking, found shortly after the laser discovery in the 
early 1960's Since then, mode- locked lasers became a 
leading way to produce ultra short pulses reaching today 
a few femtoseconds, or about two light-wave cycles. In 
a mode-locked operation, many axial modes in a broad 
frequency bandwidth are phase locked and thus provide 
one or multiple pulses in the cavity, giving at the output 
a light pulse train. The understanding of the conditions 
under which a laser operates in pulsed regime rather than 
in continuous wave regime is a question of great interest, 
both theoretical and practical. This question has been 
addressed in various studies, being referred to as "the 
second threshold" (the first one being the lasing itself) 
in the earlier years [J- [!■ 13 ■ an d recently, in the context 
of laser with a saturable absorber, as the "self starting" 
problem 0, ! ©, M M E E3 

Formation of pulses in lasers relies on the interaction 
between axial modes. Such an interaction can be pro- 
vided either by rendering the system time dependent 
(modulating) or by a suitable nonlinearity in the dynam- 
ics of the system. These two methods are commonly 
referred to as "active" and "passive" mode-locking, re- 
spectively. One type of nonlinearity known to encour- 
age pulsed operation is saturable absorption. The light 
transmissivity through a (fast) saturable absorber is an 
increasing function of the (instantaneous) input inten- 
sity. The saturable absorber destabilizes the laser op- 
eration into configurations where most of the power is 
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concentrated in short pulses. In the frequency or mode 
domain the saturable absorber induces a nonlinear four- 
wave-mixing interaction between the modes, as does the 
Kerr effect, with the difference that it is dissipative rather 
than dispersive. 

The dynamics of a laser is always subject to noise. 
Beside the usual noise sources present in every physical 
system, there is the inevitable fundamental noise of spon- 
taneous emission. This noise is inherent in lasers, since 
it always accompanies coherent amplification, on which 
lasers rely. Therefore a model of a laser that does not take 
noise properly into account risks missing key features in 
the physics of a laser system. 

The majority of laser theories treat noise as a pertur- 
bation, if at all, expecting it to manifest itself as fluctua- 
tions in the laser output. However, this approach greatly 
underestimates the effect of noise. It has been recently 
shown |14j that even very weak noise (compared to the 
intracavity power of the laser) is sufficient, for example, 
to destabilize a passively mode-locked laser, revealing a 
dramatic nonperturbative effect of the noise. 

We have recently developed 0, 0] a new approach for 
the many interacting mode system, with specific empha- 
sis on aspects of pulse formation in mode-locked lasers. 
We established an analogy between the behavior of the 
electromagnetic field (the mode system in a laser) in the 
presence of noise and equilibrium statistical mechanics, 
and applied the powerful tools of statistical mechanics to 
the problem of mode-locking. In particular it was found 
that the entropy associated with the noise is an essential 
ingredient in the theory of mode locking. This approach 
gave an inherent explanation for many experimental phe- 
nomena of mode-locked lasers, especially the existence of 
a threshold and the abruptness of formation of pulses. 
Passive mode locking was identified with a first-order 
phase transition in the model statistical mechanics sys- 
tem. Many other theoretical and experimental features 
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were found, among them, hysteresis, super-heating 
and super-cooling, successive formation of multiple pulses 
in the cavity and more. 

The results of Ref. Q were based on a mean field anal- 
ysis of the mode interaction induced by the saturable 
absorber and on numerical simulations. These indicated 
a first order phase transition when the effective temper- 
ature, i.e. the noise power, or alternatively the intracav- 
ity power, is varied. The ordered phase corresponds to 
a mode locked configuration, while the disordered phase 
to continuous wave operation. 

In this paper we show that the statistical-mechanics 
treatment of mode locking suggested in Q can be con- 
cluded with an exact solution. Using a recursion relation 
for the partition function, we are able to develop a sys- 
tematic expansion of the free energy in decreasing powers 
of the number of degrees of freedom. The first two terms 
in this expansion are calculated explicitly, allowing us 
to show that there is indeed a single first order phase 
transition associated with passive mode locking. All the 
relevant thermodynamic quantities are obtained, yielding 
a complete description of mode locking in the context of 
the present model. We also present an improved mean 
field analysis, which turns out to be exact in the thermo- 
dynamic limit. 

We show that the thermodynamics depends on a single 
parameter, the strength of the saturable absorption mul- 
tiplied by the intracavity power squared, and divided by 
the noise strength, and the transition occurs when this 
parameter crosses an explicit threshold value. 

In addition, we demonstrate the universality of the re- 
sults with respect to the gain saturation mechanism (as 
long as it is slow). It is well known that gain saturation, 
that is, the inability of the laser amplifier to amplify the 
light indefinitely, has a crucial role in understanding the 
dynamics and stab ility of lasers |l8j and in particular 
mode-locked lasers jl9j |. While in Ref. Q] gain saturation 
was incorporated by introducing a constraint of strictly 
constant intracavity power, in this paper we present a 
solution with general (slow) gain saturation mechanism. 
The relation between the two approaches is similar to 
that of the canonical and grand-canonical ensembles in 
textbook statistical mechanics, and the saturable gain 
acts in manner analogous to the chemical potential. The 
universality property follows immediately from the ther- 
modynamic equivalence of ensembles. Futhermore, for a 
given gain saturation function we are able to calculate the 
dependence of the intracavity power on the external pa- 
rameters, and in particular its jump at the mode-locking 
transition, an experimentally-measurable quantity. 

As a statistical mechanics problem, the model can be 
likened to a gas of (complex) spins, with a 4 self in- 
teraction, and a global constraint of total amplitude, 
quite similar to the constraint imposed on spins in the 
Berlin-Kac spherical model [2(j. Normally the statisti- 
cal mechanics of such systems leads to simple equiparti- 
tion. Here, however, the energy of self interaction term 
is negative, and at small enough temperature, or high 



enough power the instability stemming from the self- 
interaction drives the system into a pulsed state, where 
most of the power resides in a single spin. In Fourier 
(mode or wavenumber k) representation the model is 
equivalent to a classical complex spin chain with a special 
nonlocal interaction that drives the mode locking tran- 
sition, which in this representation is a standard order- 
ing transition. Fig. ^ shows the difference between typ- 
ical mode-locked and non-mode-locked configurations in 
Fourier space, and in real space in the context of the 
coarse grained model discussed below. 

While in this work the quartic self-interaction is at- 
tributed to saturable absorption, the model can be put 
into a broader context. The quartic interaction is of gen- 
eral interest being the lowest-order nonlinearity which is 
translation, inversion and rotation invariant and local. 
This interaction itself has been very extensively studied, 
but we are not aware of previous statistical-mechanics 
studies of its interplay with a non-local power constraint. 
In some previous well-known statistical-mechanics stud- 
ies of lasers 0|, this non-local constraint did not exist 
in the model, which may explain why the mode-locking 
noise-induced threshold behavior has not been previously 
found. Since a global limitation of power is common in 
physical systems, our model may be an important proto- 
type in nonlinear optics and elsewhere. 

II. MODELS OF PASSIVE MODE LOCKING 

A. The coarse-graining approximation 

The propagation of the wave packet in a laser cavity 
is a non-equilibrium process: A stationary condition is 
characterized by constant pumping and dissipation of en- 
ergy. Noise, due to spontaneous emission is also present. 
The temporal evolution of the wavepacket envelope ip is 
given schematically by 

d t iP(x,t) = G[ilj](x,t)+r)(x,t) , (1) 

where G is the net gain, and rj is a random driving, tf> is 
defined on the interval < x < L with periodic boundary 
conditions. 

In the context of passive mode-locking, the net gain 
functional has three essential components: G — G ga i n + 
G sa + G sp . The net pumping of energy by stimulated 
emission is modelled by 

Ggain(M) =flCP)V(z,t) , (2) 

where V — j; Jq \4>{x,t)\ 2 dx is the total power in the 
cavity. The saturable gain function g is monotonically 
decreasing with positive values for small V and negative 
for large V where various losses overcome the power sup- 
plied by the amplifier. 

A necessary ingredient for passive mode locking is sat- 
urable absorption, wherein the dissipation losses decrease 
as the power increases. Unlike the saturable amplifier, 



3 




7 1 -V % 



FIG. 1: The two "thermodynamic phases": A typical non-mode-locked configuration is shown on the left column of the figure 
and a mode-locked one on the right. The real space configurations of the coarse grained model are shown on top row, see Eqs. 
(116117ft for definitions. The Fourier modes a m , defined in Eq. (I2H at discrete wavenumbers k m , are displayed on the bottom 
row as a set of phasors, each arrow representing the complex value of a Fourier mode. 



the response of the saturable absorber is fast, so it de- 
pends on the instantaneous power \i(;\ 2 rather than on the 
total power V . We choose the specific form of saturable 
absorption 

G sa (x,t)= 7s |V>(M)|V(M) , (3) 

with 7 S > 0, valid for |-0| 2 not too large [l^. Note that 
since losses are already taken into account in G ga i„, the 
saturable absorption term describes an additional "gain" 
term which is absent when ip = 0. 

The third part models the spectral dependence of the 
amplifier that has a characteristic frequency, which we 
assume lies at the center of the wavepackct. Define the 
Fourier expansion of the wavepacket 

oo 

1>(x,t) = J2 M*)e™ . (4) 

m— — oo 

Near the resonance, the pumping efficiency falls down 
quadratically in the spectral distance. The net gain for 
modes off the center of the band is therefore reduced by 

(27T771 \ ^ 
— ) b m (t) , (5) 

where j g is positive. 

The resulting gain functional can be expressed as the 
functional derivative G(x,t) = — 5H[ip) / Sip* (x) of 

m = J dx(-^\i>(x)f+ Jg \i,'(x)\ 2 ) 

+ LU^J \ij(x,t)\ 2 dx) , (6) 



with the definition U'(V) = -g(V). 

The random part rj of Eq. will be modelled by (com- 
plex) gaussian uncorrelated noise with a power of 2LT 

(t)*(x, t)rj(x', t')) = 2TLS(x - x')S(t - t'), 

(•) stands for an ensemble average. This is appropriate 
for spontaneous emission. We arrive thus at the explicit 
equation of motion for tp 

dMx,t) = --^^+r)(x,t) , (7) 

It is well known P, that the invariant measure of 
gradient flow with additive white noise, such as Eq. Q 
is 

p[i>] = z-i e -HM/(*T) ; (8) 
where the partition function is 

Z = J[d^}[dr}e-" W,{LT) ■ (9) 

The study of passive mode locking in our model has been 
reduced to the analysis the statistical mechanics system 
described by Z. 

The model which is summarized by Eqs. H6I9|) does not 
include refractive effects such as dispersion and the Kerr 
nonlinearity (imaginary terms), which are important in 
most laser systems. Nevertheless, as pointed out in a pre- 
vious work 15], p serves as the invariant measure when 
dispersive effects are included, provided that a certain 
integrability condition holds. Ref. also studied in 
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some detail the consequences of the failure of the inte- 
grability condition to hold. The purpose of the present 
work is different: By simplifying the model even further, 
we obtain a system which captures all the essential qual- 
itative features of the passive mode locking transition in 
a transparent and intuitive manner. 

The spectral filtering term 7 ff |?/>'(a;)| 2 in H introduces 
correlations between the electric field in neighboring po- 
sitions. It counteracts the tendency of the saturable ab- 
sorption to concentrate the power in a thinner and thin- 
ner interval by introducing a length scale over which the 
electric field is correlated. This suggests a coarse graining 
approximation, in which the electric field is assumed to 
be constant on an interval, and the interaction between 
adjacent intervals is neglected. The function tj; is then 
represented on the interval by a single degree of freedom 
(see Fig. [I}. 

Taking N such intervals we find that 
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where now 
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(10) 



(11) 



(12) 



The statistical mechanics problem defined by Eqs. IlUt 
I12fl is the main object of study in this paper. At zero 
temperature, which corresponds to noiseless dynamics, 
the quartic term in H pushes all the available power into 
a single degree of freedom (whose identity depends on 
the initial conditions). This is the mode- locked state. In 
the opposite situation at high T or 7 S = the power is 
randomly distributed, and mode-locking is absent. 



B. The thermodynamic limit 

Our analysis relies crucially on N being very large. In 
the way we stated the problem, N is the ratio between 
the laser cavity length and the width of a pulse. In short- 
pulse laser this parameter is large by definition, with val- 
ues anywhere in the interval 10 2 -10 9 . Taking N — > oo 
and expanding in 1 /N looks therefore promising. 

It is not difficult to see, however, that the statistical 
mechanics system in Eqs. <|10I12H has no well defined 
thermo dynamic-limit behavior as N — > oo. For exam- 
ple, the T = value of the intracavity power V is the 
minimum of 
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~/ s NLV 2 - LU{V) 



(13) 



with respect to V, which diverges with N. As usual in 
such cases, this means that we should "renormalize" the 



parameters of the Hamiltonian, i.e. let 7 S , T and the 
function U depend on N in such a way that a well-defined 
thermodynamic limit is obtained. A comparison with a 
concrete example or experiment requires a preliminary 
translation of the renormalized parameters to the actual 
parameters. 

However, a more fundamental difficulty remains: Due 
to the nature of the mode-locking transition, the ratio 
of the peak power, the maximal value of | V'n | 2 5 to the 
intracavity power diverges like N when N is large. This 
means that only one of these quantities can achieve a 
well-defined thermodynamic limit, while the other quan- 
tity has to be rescaled. For example one can choose pa- 
rameters such that the peak power and NV achieve finite 
thermodynamic limiting values. This is convenient when 
a more general form of saturable absorber is assumed 1171 . 

In the present paper (as well as in Ref. [1|), because 
the nonlinearity is a homogeneous function (of degree 4) 
it is more convenient to choose parameters such that the 
total power V reaches a well defined limit P, and the 
peak power diverges linearly in N. When this is used in 
Eq. (|10fl the quartic term is seen to scale linearly in N. 
We choose the saturable gain function such that it scales 
in the same way, 



U N (V) = NTu(V) 



(14) 



for some fixed function u (T is added to the definition 
for later convenience). In this scheme T and j s need not 
be renormalized. 

Since the mode-locked configuration is characterized 
by peak power proportional to N, a natural order pa- 
rameter is 



1/4 



(15) 



where () stands for expectation with respect to the in- 
variant measure. The peak power in the non-mode-locked 
case is 0(1) in the thermodynamic limit, so M = there. 



C. Fixed power ensemble 

Since the quartic term in the Hamiltonian H is un- 
bounded from below, the gain saturation term U (V) is 
essential to ensure stability, preventing the system from 
cascading into states with arbitrarily low H . This re- 
flects the well known fact that lasers owe their stability 
to gain saturation [181 [l9j . This is analogous to the role 
of the chemical potential in the grand-canonical ensem- 
ble in standard statistical mechanics, which limits the 
number of particles in the system. 

An alternative approach is to suppose that the intra- 
cavity power V has a fixed value P, whose analogue in 
textbook statistical mechanics is the canonical ensemble 
where the number of particles is fixed. This is the scheme 
used in Ref. One of the main results of the present 



5 



work is an equivalence of ensembles. The thermodynam- 
ics obtained in the fixed- and variable-power ensembles 
are equivalent. The variable-power ensemble has to be 
used if (V) is not known. These issues are discussed in 
sec. E| 

In the fixed power ensemble there is no need to include 
the gain saturation term, and the partition function is 
defined by 



with the reduced Hamiltonian 
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N 



(16) 



(17) 



The change of variables y n = \ip n \ 2 in Eq. (|16f) leads to 
the simpler form 

Z N ( ls ,T,P)= l'l[dy n e^Z n y 2 n5 (l_'£y n _p\ , 

^ n \ n=l / 

(18) 

where the y n integrations are from to oo. Another 
change of variables y n — > Py n leads to the useful scaling 
relation 



Zjv( 7s ,T, P) = P N - 1 Z N ( ls P I /T) 



where 



(19) 



An important conclusion has already been reached: The 
thermodynamics depends on the single parameter 7 = 
7 S P 2 /T. Eq. (|19fl proves this in the fixed power scheme, 
while the equivalence of ensembles extends this to the 
general case. 

In Sec. IIVI we solve the statistical mechanics problem 
in the fixed-power ensemble by developing an asymptotic 
expansion for Zm for large N, which is valid and uniform 
for all 7. 



III. MEAN FIELD THEORY 

The problem lends itself to a mean field analysis when 
formulated in Fourier (mode) space. In the Fourier rep- 
resentation mode-locking manifests as ordering in the 
phases of the various modes, see Fig The discrete 
Fourier transform of the Hamiltonian (|17fl is 



mi —m^-\-m-i—va^—^)N 



where 



N 



4>n = J2 a™ e2 "" lt 



(20) 



(21) 



and p is an integer (actually it is easy to see that only 
p = — 1, 0, 1 are relevant). V is now given by 



v = E k 



(22) 



The main disadvantage of the Fourier space formulation 
is that the nonlinear term becomes complicated and non- 
local. Mean field theory overcomes this difficulty by as- 
suming that the different modes are uncorrelated and 
characterized by a common probability distribution func- 
tion /9 m f(a). When the problem is formulated in Fourier 
space, the mean field approximation is quite plausible, 
since the interaction term involves all the degrees of free- 
dom. 

In the mean field framework the free energy per degree 
of freedom is 



F 



\ogZ_ 

N 



2T 



N 2 \(a} ml \ 4 + (logp) miJ (23) 



where () mf stands for expectation value with respect to 
p m f. Gain saturation is included by demanding that 



<M 2 ) mf 



p 

N 



(24) 



Following the standard procedure of mean field calcu- 
lations [23J, p m f is found by minimizing the free energy 
subject to the constraint Ea. l24|) . A necessary condition 
for the minimization of F is stationarity with respect to 
variations of p, 



= 



Sp(a] 



(F + \((\af)-P/N)) = 



(25) 



= -2Re(-f s /T)N 2 \(a)\ 2 (a*)a+ (logp(o) + 1) + A|a| 2 , 

where A is a Lagrange multiplier. The solution of Eq. 125|) 
is a gaussian probability distribution function 



1 |q-<q)| 



P( a ) = 2 e 

ira 

(a) and a are related by Eq. I|24J) which implies 



<I« |2X 



/mf 



ml' 



P 

N 



(26) 



(27) 



p m { is therefore characterized by the single parameter 
M = \/N/P\{a) m{ \. From Eq. {23 one can see that 
< |M| 2 < 1. When M = the phases of the modes are 
completely random, which means that in real space the 
power is uniformly distributed. When M > the modes 
are correlated which in real space means that a macro- 
scopic fraction of the power resides in the variable i/jn- 
M is therefore an order parameter, which can be shown 
to asymptotically coincide with the previous definition of 
the order parameter Eq. (|15J) . 

We note that the present formulation of mean field 
theory is slightly careless, in that it misses ordered states 
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centered at other points than i/j^ in real space. The in- 
clusion of these configurations would require considering 
cases where various a m 's differ by a fixed phase. How- 
ever, since the number of missed configurations grows 
only linearly with N, the entropy is underestimated by 
a subextensive factor, and the thermodynamics is not 
affected. 

We can define a new thermodynamic potential /(7, y), 
where y = M 2 and 7 = j s P 2 /T, which is the free energy 
for a given value M. Using Eqs. I|26l27jl in Eq. (J2HJ it is 
found that 

f(i,y) = -(ly 2 + iog(i-y)) , (28) 

up to unimportant additive terms independent of 7 and 
M. The free energy F for given 7 is the global minimum 
of /(7, •), and the abscissa of the minimum, y, is the 
square of the order parameter. 

The function / of Eq. I|28[l has a single minimum for 
7 < 4, at yo = 0. The vanishing of M means that 
the phases are not locked corresponding to the disor- 
dered, non mode-locked phase. For 7 > 4 there exists 
an additional (local) minimum, y\ > 0, see Fig. IIIII 
which corresponds to the mode-locked state. However, 
for 4 < 7 < 7* ~ 4.91 f{j,yi) > f(j,yo), which means 
that the mode-locked state is metastable, and the true 
equilibrium is still disordered, y = yo. At 7 = 7* the 
two minima exchange stability, and for all 7 > 7* the 
equilibrium state is mode locked, with ^fy{ as the value 
of the order parameter. 7* is the solution of the equation 
/(T.J/o) = /(7> 2/1(7))) 




FIG. 2: The thermodynamic potential f(-y, y) as a function 
of y = M 2 for several values of 7 ranging from 2 to 6 with 
higher values of 7 corresponding to lower curves. Curves with 
two (local) minima correspond to systems with a metastable 
state. The critical 7 (7*) is by definition the one where the 
values of f(-y, ■) at the two minima are equal, which implies a 
first order phase transition. 

In terms of the original variables 7 S , T and P, the 



phase transition point is therefore 

\~ = 7* « 4.91 (30) 

In order to compare it to the result in Ref. , one should 
remember the difference in the modeling of spectral fil- 
tering here and there. In Ref. [l| we imposed "Dirichlet" 
boundary conditions in the Fourier space, while here, the 
coarse graining method actually induces periodic bound- 
ary consitions in Fourier space. Since the interaction in 
Fourier space is long ranged, this leads to a difference. 
The Hamiltonian in Ref. [lj was identical to Eq. I|20|) , ex- 
cept that p — only. The number of quartets with k = 
is 2/3 of their number with p = — 1,0, 1. In the mean 
field approximation this would simply lead to a factor of 
3/2 in the transition temperature: 

-v P 2 3 

^- = ^7*«7.4 (31) 

This is very close to the result of j s P 2 /T w 7.7 obtained 
in Ref. 0. (We note that the mode locking transition 
was specified in Ref. pj in terms of I/7.) The mean field 
theory presented here is better, since here we do not make 
the ansatz of separating the modulus and angle in p m { . 
Here we also end up with the simple analytical expression 
H28JI . Since the Hamiltonian l|20Jl is translation invariant 
in Fourier space, mean field theory is more suitable for 
its analysis than for that of its counterpart from Ref. . 

In summary, it has been demonstrated in the mean 
field context that the mode locking transition is a stan- 
dard first order transition, accompanied by coexistence 
and metastable configurations in its neighborhood. The 
mean field theory involves an uncontrolled approxima- 
tion, which is hard to justify rigorously. In the present 
work the justification will ultimately follow from the real 
space analysis given in the next section. 

Evidently, the transition condition depends, although 
not greatly, on the spectral filtering scheme. Analysis of 
the parabolic filtering scheme, which we do not pursue 
here, yields a yet another value for the transition tem- 
perature, close to (|3TJ)l and il3*T|) . 



IV. FIXED POWER FINITE N ANALYSIS 

In this section we calculate an asymptotic expansion 
of the partition function Z/v(7, P) in decreasing powers 
of N, the number of degrees of freedom. All information 
pertaining to the passive mode-locking transition is then 
obtainable in a standard manner. In particular, we show 
that the mean field calculations give the exact free energy. 

Our starting point is a recursive version of Eq. i|18fl for 
Zjv, obtained by performing only N — 1 of the y integra- 
tions, 
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„ 2 / N — 1 iV — Vat 

Z N ( ls ,T,P)= I dy N e*frrV"Zs-i[^-7.,T, - m P 



N 



N - 1 



(32) 



After using the scaling relation Eq. (|19|) and making a change of the integration variable we obtain a recursive equation 
for Z: 



(33) 



This is the fundamental equation of the real space anal- 
ysis. 

We will show that when TV is large the only significant 
contribution to the y integration in Eq. I|33() comes from 
the vicinity of one or two values of y which maximize 
the integrand, one of which is yo — 0. The integration 
on other parts of the interval is exponentially small in N 
and will be neglected. The case of a single maximizing 
point will be shown to correspond to invariant measures 
concentrated on configurations where the amplitude of 
all degrees of freedom is 0(1), i.e., non-mode locked, dis- 
ordered configurations. This happens for small enough 7. 
When there are two maximizing points the typical con- 
figurations are such that a finite fraction of the power is 
concentrated in a single degree of freedom, while the am- 
plitude of other degrees of freedom is again O(l). These 
mode locked configurations arise for large enough values 
of 7. We show that these are the only two possibilities. 



A. The disordered phase 

We first tackle the case of small 7. To this end we 
use the Fourier representation of the delta function in 
Eq. i fTSjl to reexpress Z by 



Zn{h) 




N 



TV +Jj z V 



(34) 



for some Cn > N. As long as Cat > A, the right- 
hand-side of Eq. (|34|l is independent of Cn- One can 
now expand the quadratic term in the exponential in a 
Taylor series keeping the first two terms, carry out the y 
integration, and then take Cn — ► 00 giving 



Zn{i) 



J 2ni V z 



N 



(35) 



The contour of integration has to be deformed so as 
to avoid the singularity at z — 0. A standard argument 
shows that the contour should be moved to the left so 
that it crosses the real line at a negative value. The inte- 
gral (|35|l can then be calculated by pushing the contour 
through the singularity at z = and then to Rez = 00. 
The exponential in the integrand makes the integration 
at infinity vanish in the limit, leaving only integration on 



a contour surrounding z — clockwise. Using Cauchy's 
theorem this evaluates to 



Zn{i) 



N N 

E 



dz ^/ N \ (jNT 



n\(N - n)\(N + 2n)\ ' 
or, using Stirling's formula, 



(36) 



Zn{i) 



zffi certainly provides an asymptotic approximation 
of Zn{^) as 7 — ► 0, but we shall show that it also serves 
as the leading term of Zn^) as N — > 00 for all < 
7 < 7*. This is achieved by showing that the recursive 
equation 

4 0) (7)~^ / dyel^d-yf^Z^d-yf) 



= 4 0) ( 7 ). (37) 



= N- 



V2nN Jo 



(1 - y) 



(38) 



holds for 7 in this interval. Here and below the sym- 
bol ~ stands for 'asymptotic for large N\ Consider 
the last integral. When N — > 00 the integrand be- 
comes strongly peaked near the global minimum 2/(7) of 
f(l>y) = — ( 2y 2 + l°S(l — y))- The function / is precisely 
the thermodynamic potential encountered in the context 
of the mean field approximation, see Eq. (J2SJ. As shown 
above fSec. lIIIf) . 3/(7) = yo = for 7 < 7*. For these val- 
ues of 7 only the neighborhood of y = has to be taken 
into account and the integral in the right hand side of 
Eq. I|38|) is approximately 



Ae 7 - 



V2nN Jo 
which establishes 



dy (1 - y) 



N-2 „ e 7 



^(7) ~4 0) (7) , 7<7* 



Z ( n\ (39) 



(40) 



The thermodynamics now follows straightforwardly. 
For example, the free energy per degree of freedom is 



8 



F(-f) = -^logZAr ~ 1, independent of 7 to leading or- 
der, and the expectation value of |^| 4 in the invariant 
measure is 



4\ _ 



(7 < 7*) , 



(41) 



also independent of 7 in the leading order. In particular, 
the order parameter M from Eg. 1151 equals zero, showing 
that this is indeed a disordered configuration. 



B. The mode-locked phase 



We turn now to the case 7 > 7*, where y = y\ > 0. 
We can no longer expect that Zn ~ Zj!?\ but the mean 



field calculations suggest that 



Z N {l)~ A N {l)e~ NF ^ 



(42) 



where ^(7) = / (7, 2/(7)) — 1 and An is subexponential 
in N . The results of the previous subsection imply that 
the asymptotic form 14211 is valid for 7 < 7*, since then 
y = yo, and F = —1. Using Eq. l(3*3l for 7 > 7*, we 
presently show that Eq. Ij42(l is valid for all 7 7^ 7* and 
find explicit expressions for An- 

Substituting Eq. 142|) in Eq. i.YMi gives the asymptotic 
equation 



J 



i-N 



( 7 ) e -JVF( 7 ) r~ Ne dy j 



e 2 1 



-.4 



JV-l 



(7(1 - y) 



2\ e ~(N-l)(f^ : y)+F( 7 (l-y) 2 )) 



(43) 



r 



As before, the integration is concentrated near the max- 
imal points of the large exponential, i.e. the minima (as 
a function of y) of f(j, y) + F(j(l ~ J/) 2 ))- Recalling the 
definition of F the minimization problem turns into 

min/( 7 ,y)+F(7(l- 2/ 2 )) = 
y 

= min— 7y 2 — log(l — y) — 7(1 — y) 2 w 2 — log(l — w) = 

y,w 



7 

min - o (y 2 + w2 ) + ^sC 1 - v - w )' 



(44) 



y,w 2 

where we have put w — w(l — y). It is straightforward 



to check that either j/orro must vanish at the minimum. 
For 7 > 7* there are two possibilities, y — 2/1(7), w = 0, 
and y = 0, w = 2/1(7), and the minimal value is the same 
in both cases. 



The conclusion is that the integration receives two 
main contributions, one from the neighborhood of y = 0, 
and one from the neighborhood of y = y, which we denote 
by I and I\ respectively. Iq is found by approximating 
the exponential near y = and evaluating prefactors at 
y = giving 



J 



Jo~iV e A w _i(7)e-( JV - 1 ) F W / dy{\-y) N 

Jo 



-2 e 2(N~l)F'( 7 )y 



1-2^(7) ™ 



(45) 



r 



the assumption that An is subexponential in N was used 
to approximate An-i ~ An- 

For the calculation of I\ we need to evaluate F near 
7(1 — y( 7 )) 2 . It follows from the properties of / that this 
is always strictly less than 7*, where F = — 1. Therefore 



,NF(n) 



2tt(1-2/) 2 

e 7(l-S)' 



dye 



(N-l)F^>(j) 



(2) ,..\ (y-vr 



-NF(~f) 



(i-y) 2 VF^h) 



(46) 



where F^'{"f) = <9 2 /( 7 , 2/(7)) ■ Combining these results 



with Eq. I|43(l . we get a linear equation for An whose 
solution is 



1 - 2F'( 7 ) 



(l-y) 2 V / F( 2 Tl-2F'(7)-e^)+ 



l 7 > 7 



(47) 

This establishes equation Eq. I|42|) , with explicit values 
for A N for all 7^7*. 

We find now that -F(7) is indeed the free energy, con- 
sistently with the mean field theory. Since 



<|^| 4 )~2(logZ w )'(7)-2 



^(7) 
A( 7 ) 



2iVF'( 7 ) , (48) 
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the order parameter M = (2^(7)) 1 4 is nonzero for 7 > 
7*, showing that mode locking occurs for such 7. Using 
the definition of F we can calculate explicitly 

-2F'( 7 ) = ~f{i,yin)) = -d-yf{i,y(i)) = f , 

7 . (49) 

since / is by definition stationary with respect to y at y. 
This means that M = y/y, also in accordance with the 
mean field calculations. It is possible to show by calculat- 
ing higher moments that M 2 is the power concentrated in 
a single degree of freedom. This result, which in physical 
terms means that mode locking results in a single pulse, 
can be traced to the fact that the integral in Eq. 143(1 
receives contributions only from yo and y\. 



C. The transition region 

The analysis of the previous section does not apply to 
the case where 7 is precisely equal to 7*. For example, 
Eq. 1471) would imply that Ajv(7*) is infinite, since the 
derivative in F'{^) should be taken from below. More 
importantly, the asymptotic approximation Eq. 1(42(1 is 
not uniform in N near 7* , because it neglects the contri- 
bution of the metastable state. An asymptotic approxi- 
mation for Z which is valid and uniform for all 7 is 



N 



+ A* N {l)e 



-NF'(j) 



(50) 



where F*(j) = f(j, 2/1(7)), and A* N is given by Eq. lET7|) 
replacing everywhere F by F* and y by y±. The uni- 
form approximation is a continuous function of 7 which 
reduces to the non-uniform approximation for (7 — 7*| 3> 
1/N. 

Observables in systems with a finite number of degrees 
of freedom exhibit crossover behavior in the mode- locking 
transition, rather than the sharp, discontinuous depen- 
dence on parameters predicted in the thermodynamic 
limit. When the number of degrees of freedom is not 
too large, the crossover is measurable and describable by 
the uniform approximation. For example, 



(\tP\ i )~2(a N ( 1 )+Nb N ( 1 )y 1 ( 1 ) 2 ) 



where 



and 



7 T7&F + (^)'(7) e - A ' F ^ ) 



a N 



J^L 1 A* , 



-NF'(j) 



b N = 



^(7)e- JVF ' (7) 
£Lr + A* M e~ NF *(~f) 



(51) 



(52) 



(53) 



For moderate values of N, there is a significant interval 
in 7 below the transition where (IV'I 4 ) is much larger 
than its value in the thermodynamic limit. A comparison 
between the uniform and non-uniform approximations to 
(\t(j\ 4 ) is shown in Fig. llVH 




4.4 4.6 4.8 5 5.2 5.4 



Y 



FIG. 3: A comparison between the value of the order param- 
eter in the uniform (full line) and the non-uniform (broken 
line) approximation as a function of the nonlinearity param- 
eter 7 for a system with N = 150 degrees of freedom. The 
non-uniform approximation breaks down near the transition 
point 7*. 



V. THERMODYNAMICS WITH VARIABLE 
TOTAL POWER 

In the previous section we developed a systematic ap- 
proximation scheme for the partition function Z of the 
passive mode locking model as a function of the nonlin- 
earity strength j s , the fixed total intracavity power P 
and the temperature T. This allowed us to calculate the 
free energy energy per degree of freedom, and we found 
that mode locking occurs whenever r y s P 2 /T is greater 
than a critical value 7*. 

However, in experimental situations the intracavity 
power V is not fixed in advance. Rather it is a fluctuat- 
ing quantity, whose mean value P is determined by the 
saturable gain function U (see Eq. iJTJJl). The relation be- 
tween the thermodynamics in the fixed-power ensemble 
analyzed above, and the variable-power ensemble which 
is the subject of this section is quite similar to the one 
between the canonical and grand canonical ensembles in 
statistical mechanics j2^. In the latter case one defines 
the grand potential £1 = fiN — F, where /i = dF/dN is 
the chemical potential. An equivalent thermodynamics 
is obtained after replacing the extensive variable N, by 
the intensive variable \x. Finite size corrections to the 
thermodynamics in the two ensembles are also related, 
but not equivalent. 

In this section we show in a similar spirit that ther- 
modynamics with fixed and variable power is equivalent, 
and calculate the thermodynamic limit of P = (V) in the 
variable power case. It is quite straightforward to gener- 
alize the calculations and to obtain subleading terms as 
in Sec. IIVI but this is not pursued here. 

In the fixed power ensemble the free energy per degree 
of freedom F of the Hamiltonian Hjy which includes the 
saturable gain is related to the free energy F calculated 
in Sec. E] by 



F(j s , T, P) = F{ ls P 2 /T) + log P + u(P) 



(54) 
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(Refer to Eas. (|10ll4(l for the relevant definitions). We 
now define the variable-power thermodynamic potential 
$ as the Legendre transform of F, 

$( 7s ,T, fj) = mm/xP - H 7s ,T, P) . (55) 

The function u has to grow faster than P 2 to ensure that 
the minimum in Eq. (|55|l is finite. We also require that 
u is convex to obtain a unique minimum. Otherwise u is 
arbitrary. 

The thermodynamics is obtained from $ through the 
properties of the Legendre transform, the case of interest 
being fi = 0. The mean power, 



P = d^( 7s ,T,0) 



(56) 



can be found from the definition Eq. (|55|l and the results 
of the previous section; it is given implicitly by 

1 + Pu'(P) - 7 y( 7 ) 2 = , (57) 

where as before 7 = jsP 2 /T. The order parameter is 

M 4 = 2T<9 7 $ = -2Td y F = -P 2 F'{ 1 ) = P 2 y(-/) 2 , 

(58) 

which, for a given mean power P, is independent of the 
form of g, and therefore also equal to the order parame- 
ter in the fixed power ensemble. Moreover, the thermo- 
dynamics depends on the single parameter 7. A special 



case is the mode-locking transition point, which occurs at 
7 = 7* whatever the form of the saturable gain function. 
This universal behavior stems from the thermodynamic 
equivalence of the fixed power and variable power ensem- 
bles. 

Another interesting thermodynamic quantity which 
can be studied only in the variable power framework is 
the susceptibility \ = P'( 7 ) which measures the response 
of the intracavity power to changes in the strength of the 
nonlinearity or inverse noise power. Taking the derivative 
of Eq. JS7J shows that 



(Pu'(p)y 



(59) 



In the non-mode locked regime P is independent of 7 and 
X = 0; when y > and mode-locking occurs, j/'( 7 ) is also 
positive and it follows from the convexity of u that the 
susceptibility is strictly positive in mode-locked systems. 
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